Transport and bistable kinetics of a Brownian particle in a nonequilibrium 

environment 
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A system reservoir model, where the associated reservoir is modulated by an external colored 
random force, is proposed to study the transport of an overdamped Brownian particle in a periodic 
potential. We then derive the analytical expression for the average velocity, mobility, and diffusion 
rate. The bistable kinetics and escape rate from a metastable state in the overdamped region are 
studied consequently. By numerical simulation we then demonstrate that our analytical escape rate 
is in good agreement with that of numerical result. 
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I. INTRODUCTION 

As an immediate consequence of stochastic dynamics, 
it is observed that thermal diffusion in a periodic po- 
tential plays a prominent role in various cases such as 
Josephson's junction^, diffusion in crystal surface^, noise 
limit cycle oscillators^ etc. There has been a renewed 
interest in recent times in the study of transport proper- 
ties of Brownian particles moving in a periodic potential^ 
with special emphasis on coherent transport and giant 
diffusion^. These studies have been motivated in part 
by an attempt to understand the mechanism of move- 
ment of protein motors in biological systemsS.. Several 
physical models have been proposed to understand the 
transport phenomena in such systems such as vibrational 
ratchet-, rocking ratchet^, diffusion ratchet^do^ correla- 
tion ratchetii, etc. Such ratchet models have a wide 
range of application in biology and nanoscopic systemsi^, 
because of their extraordinary success in exploring exper- 
imental observations on biochemical molecular motors, 
active in muscle contractioni2., observation of directed 
transport in photovoltaic and photoreflective materials^^, 
etc. In all the above models, the potential is taken to 
be asymmetric in space. One can also obtain a unidi- 
rectional current in the presence of spatially symmet- 
ric potential. For such non-equilibrium systems, one re- 
quires time asymmetric random force- - or space depen- 
dent diffusio n^^d^d^ . In passing, we want to mention 
the fact that to explain the role of Levy stable noise in 
transport phenomena in presence of bistable, metastable, 
and periodic potentials with broken symmetry, several el- 
egant approaches have been suggested recently 

Traditionally, Langevin equation describing the dy- 
namics of a Brownian particle coupled to a thermal reser- 
voir is a tool for modeling several aspects of nonequilib- 
rium phenomena^. In addition to the Langevin dynam- 
ics, one often takes into account the probabilistic aspect 
of the random dynamics through the usage of Fokker- 
Planck equation^^. Both of the approaches utilize the re- 



lation between the random fluctuations imposed by the 
reservoir into the Brownian particle and the relaxation 
of the imposed energy back to the reservoir, through the 
fluctuation-dissipation relation (FDR)^. FDR takes into 
account the balance between the energy input (to the sys- 
tem from the reservoir) and output (from the system into 
the reservoir) through the detailed balance mechanism. 
Typical signature of such a system is the attainment of 
equilibrium in the asymptotic limit. 

An additional external random driving applied to the 
Brownian particle can break this balance mechanism and 
make the composite system open-^, the direct conse- 
quence of which is the loss of FDR. In addition to that, 
the system hardly reaches the equilibrium state in the 
long time limit but rather attains a stationary steady 
stated?. However, a driving of the reservoir by an exter- 
nal random field^i creates a thermodynamic consistency 
condition analogous to FDR3^ that leads to the study 
of several interesting phenomena3^iSSii^'22ii2S in chemical 
physics. The net effect of the reservoir driving by an 
external random force is the creation of an effective tem- 
perature in addition to the thermal energy fc^T exerted 
by the reservoir on the system of interest. As shown 
recently, this effective temperature can enhance the re- 
action rate in condensed medicfc^^i^SiiSl as well as generate 
directed motion in a periodic syste m^^'^^ . 

In the present paper, we consider a system-reservoir 
model where the bath is modulated by an external noise. 
However, when the reservoir is modulated by an external 
noise it is likely that it induces fluctuations in the po- 
larization of the reservoir due to the external noise from 
a microscopic point of view and one may expect that 
the non-equilibrium situation created by modulating the 
bath makes its presence felt in the transport property and 
also in the kinetics of the Brownian particle. A number of 
different situations depicting the modulation of the bath 
may be physically relevant. For an example, we may con- 
sider the case of a Brownian particle when the response 
of the solvent is time-dependent, as in a liquid crystal, 
or in the reaction-diffusion mechanism in super-critical 
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lattice, or the growth in hving polymerization^^. 

To observe the effects of external stochastic modula- 
tion one can carry out the experiment in a photochemi- 
cally active solvent (the heat bath) where the solvent is 
under the influence of external monochromatic light with 
fluctuating intensity which is absorbed solely by the sol- 
vent molecules. As a result of this, the modulated solvent 
heats up due to the conversion of light energy into heat 
energy by radiation less relaxation process, and an effec- 
tive temperature-like quantity develops due to constant 
input of energy. Since the fluctuations in light intensity 
result in the polarization of the solvent molecules the ef- 
fective reaction field around the reactant gets modified'^^. 
Our theoretical model can be tested experimentally to 
study the directional motion and mean first passage time 
of artificial chemical rotors in photovoltaic solvent^. 

There are some precedents for our model that are 
worth mentioning. Bravo et al^ dealt with a related 
problem: a classical system in a heat bath with an addi- 
tive external noise. In the quantum case, Faid and Fo5s^ 
proposed a stochastic coupling between the system and a 
heat bath as a phenomenological mechanism of relaxation 
for the bath. Manas et al.^"^ considered the problem of a 
system coupled to an ensemble of independent harmonic 
oscillators as a reservoir. We have also considered the 
dynamics of a metastable state non-linearly coupled to 
a heat bath driven by an external noise to study the es- 
cape rate from a metastable state^^. In the present work, 
we address the so-called ratchet problem and bistable ki- 
netics of a Brownian particle for a thermodynamically 
open system where the associated bath is modulated by 
a colored noise and we explore the dependence of var- 
ious parameters of the external noise on the transport 
phenomena and bistable kinetics. 

To get insight of situations encountered in the grow- 
ing number of nano and microscale experiments, our 
model can be used as a potential tool. For instance, sev- 
eral DNA nanomotors have been recently suggested^^*^. 
These machines are relatively slow and do not per- 
form continuous rotation. Very recently, a rotary DNA 
nanomachine that shows a continuous rotation has been 
proposed—. This motor consists of a DNA ring whose 
elastic features are tuned such that it can be externally 
driven by a periodic temperature change. Our model pro- 
posed in this paper can be used as a theoretical avenue 
to examine the periodic temperature change via a phys- 
ically motivated microscopic Hamiltonian picture. Opti- 
cal tweezers^l, which is capable of manipulating nanome- 
ter and micrometer sized dielectric particles by exerting 
extremely small forces via a highly focused laser beam, is 
often used to manipulate and study single molecules by 
interacting with a bead that has been attached to that 
molecule. DNA and the proteins and enzymes that in- 
teract with it are commonly studied in this way. At this 
juncture, we want to mention the fact that our develop- 
ment present in this article can be used as a theoretical 
model to understand the phenomena of heating of the 
liquid surrounding a bead in an optical tweezers setup^. 



The organization of the paper is as follows: In Sec. 
im starting from a microscopic Hamiltonian picture of a 
system linearly coupled with a harmonic reservoir which 
is modulated by a noise with arbitrary decaying memory 
kernel, we have derived the Langevin equation with an 
effective noise ri{t) and then explored its statistical prop- 
erty. Employing the functional calculus method2^>2Si^, 
we then obtain the Fokker-Planck-Smoluchowski equa- 
tion in Sec. Illli corresponding to the Langevin equation 
valid in the over damped limit and for rapid fluctuations 
whose correlation function vanishes rapidly. In Sec. IIVI 
we have calculated the steady current in a ratchet po- 
tential and derived the expression for diffusion rate and 
mobility. As another application of our development, we 
study the bistable kinetics to obtain the stationary prob- 
ability density function (PDF) and the barrier crossing 
rate. The summarizing remarks are presented in Sec. IVII 
preceded by a numerical application in Sec. [V] 



II. THE MODEL: HEAT BATH MODULATED 
BY EXTERNAL NOISE 

We consider a classical particle of unit mass coupled 
to a heat bath consisting of a set of A'^-numbers of mass 
weighted harmonic oscillators with frequency {ujj}. The 
heat-bath is externally driven by a Gaussian random 
force e{t) with an arbitrary decaying correlation func- 
tion. The total Hamiltonian of such a composite system 
can be written as^^''''^ 



H 



V{x) 



N 

E 



Cjxf 



i^int.(l) 



In the above equation x and v are the coordinate and the 
velocity of the system particle, respectively, and V{x) is 
the potential energy of the system, {xj, Vj} are the vari- 
ables for the jth oscillator with characteristic frequency 
LDj. The system particle is coupled to the bath oscillator 
linearly through the general coupling terms CjtOjX where 
Cj is the coupling strength for the system-bath interac- 
tion. The interaction between the heat bath and the 
external noise is represented by the term Hi^t which we 
take afM^ 



N 

Hint = y^^KjXjejt), 



(2) 



where Hj denotes the coupling strength of interaction 
and e(t) is an external noise which is assumed to be sta- 
tionary, Gaussian with zero mean and arbitrary decaying 
correlation function, the statistical property of which is 
given by 



(e(i))e = 0,(e(t)e(t'))e -2Z?,*(i-t'), 



(3) 



where is the external noise strength and '^{t) is the 
external noise memory kernel which is assumed to be 
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a decaying function of its argument and (■ • • )e implies 
averaging over each realization of e{t). Eliminating the 
bath degrees of freedom in the usual way^^^, we get the 
Langevin equation for the system particle 



V'{x) - / dt'jit - t')v{t') + f{t) + iT(t), (4) 



where the memory kernel 7(t) and the Langevin force 
term f{t) are given, respectively, by 



N 



f{t) = ^c^2w|cos(a;j0, 



(5) 



i=i 

N 



sin(a;_,i) 



(6) 



In Eq.(j4]), Tr{t) is the fluctuating force generated due to 
the external stochastic forcing of the bath by e{t) and is 
given by 



At) 



dt'ip{t - t')e{t'), 



with 



N 



ip{t) — CjLUjKj sm{ujjt). 



(7) 



(8) 



The form of Eq.(j4|) reveals that the system is driven 
by two fluctuating forces, f{t) and 7r(t). n{t) is a dressed 
noise originating due to the bath modulation by external 
noise e{t), and f{t) is the thermal noise due to system- 
bath coupling. To deflne the statistical properties of f{t), 
we assume that the initial distribution is such that the 
bath is equilibrated at t = in the presence of the system 
but in the absence of the external noise e(t) such that 



(/W> =0, {f{t)f{t'))^kBTj{t-t'), 



(9) 



where ks is the Boltzmann constant and T is the equi- 
librium temperature, (• • • ) implies the usual average over 
the initial distribution which is assumed to be a canonical 
distribution of Gaussian formS^i^ 



P = 7Vexp' 



2kBT 



where Af is the normalization constant. Now at t = 0+, 
the external noise agency is switched on to modulate 
the bath. Here, we define an effective Gaussian noise 
r]{t) = f{t) + 7r(t) the statistical property of which can 



be described by 

{{vim 

{{vitMt'))) 



= 0, 



dt" 



ksT-fit - t') + 2De I dt" I 
Jo Jo 

Xip(t - t")ip{t' ~ t"')^{t" - t'"), (10) 
G{t - t') (say). 



In Eq. fTOl) . ((•••)) means that we have taken two av- 
erages independently, average over initial distribution of 
bath variables and average over each realization of e(t). 
While deriving Eq.dTU]), we have made the assumption 
{{ri{t)ri{t'))) = G{t — t'), which can not be proved un- 
less the structure of Lp{t) is explicitly given. As we shall 
see later it is a valid assumption for a particular choice 
of coupling coefficients and for external stationary noise 
processes. It should be realized that Eq.([TU]) is not a 
fluctuation-dissipation relation (FDR) due to the appear- 
ance of the external noise intensity, rather it serves as a 
thermodynamic consistency relation. 

To obtain a finite result in the continuum limit, i.e., 
for A'^ — > CXI, the coupling function Ci = c{u)) and Hi = 
k{ijj) are chosen as c{lo) = cq/uj^/t^ and k(uj) = kqoj^Jt^. 
Consequently, ^[t) and {p{t) reduce to 



lit) 



r2 



dojpiijj) cos wt. 



and 



</3(t) = CqKo / dLOpiljj)Ld smLot, 



(11) 



(12) 



where cq and kq are constants and is the correla- 
tion time of the heat bath. For Tc ^ we obtain a 
(5-correlated noise process. I/tc may be characterized as 
the cut off frequency of the bath oscillators. p(a;) is the 
density of modes of the heat bath which is assumed to 
be Lorentzian 



p(cj) 



2t, 



7r(l -t- w2r2) 



(13) 



The above assumption resembles broadly the behav- 
ior of the hydrodynamical modes^^''*'* in a microscopic 
system and is frequently used by the chemical physics 
communityiS. With these forms of pioj), c(u;) and k(w), 
we have the expression for ip{t) and 7(i), respectively, as 



CqKo 



exp(-|t|/r,),7(t) = ^exp(-|i|/T,). 



At) 



Though Eq.Hni) is not a FDR, Eq.® resembles the fa- 
miliar linear relation between the polarization and the 
external field. Here, 7r(i) and e(t) play the role of former 
and latter, respectively. Thus ^pit) can be interpreted as 
a response function of the reservoir due to external noise 
e(i). It is also clear from the structure of ip{t) and 7(i) 
that 



dlit) 
dt 



Co J_ 



Vit)- 



(14) 
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The above relation is independent of 'Diuj) and represents 
how the dissipative kernel j{t) depends on the response 
function of the medium due to the external noise e{t). 
Such an equation for the open system can be anticipated 
in view of the fact that both the dissipation and response 
functions crucially depend on the properties of the reser- 
voir. If we assume that e is a (5-correlated noise, i.e., 
{e{t)e{t'))e = 2De6{t — t') then the correlation function 
of 7r(i) is given by 

{At)At'))e = exp(-|(i - i')l/rc), (15) 

where we have neglected the transient terms {t,t' > Tc). 
This equation shows how the heat bath dresses the exter- 
nal noise. Though the external noise is a (5-correlated one, 
the system encounters it as an exponentially correlated 
noise with the same correlation time of the internal noise 
but with a strength dependent on the coupling term kq 
and the external noise strength De- On the other hand, 
if the external noise follows Ornstein-Uhlenbeck process 

(e(0e(O>e = — exp(-|(i-t')|/re), 

Tc 

the correlation function of n(t) is found to be 



-exp(-|(t-t')|/T,) 



(16) 



where we have again neglected the transient terms. The 
dressed external noise 7r(i) now has a more complicated 
structure of correlation function with two correlation 
times Tc and Tg. If the external noise correlation time 
is much larger than the internal noise correlation time ( 
Te » Tc), which is more realistic, the dressed noise is 
dominated by the external noise and we have 



that 

-fit) = 2cl6(t - t') = 2-fS{t - t') 

(fit)) = 0, 

{fit) fit')) - 2clkBTS{t-t') = 2jkBT5it^t'), 

with 7 = Co [see Eqs.(l5|), (O and ([7|)]. The effective noise 
riit) thus has statistical properties [see Eq.®] 



Mt))) 

imvit'))) 



0, 

Dr 

tr 



exp(-|(t-t')l/Tfl), (20) 



where 



Dr = 'fiksT + DcKq) and Tr = ^JuItc, (21) 

with Dr and tr being the strength and correlation time 
of the effective noise r/it), respectively. Although the 
reservoir is driven by the colored noise e(t) with noise 
strength Dg and correlation time Tc , the dynamics of the 
system of interest is effectively governed by the scaled 
colored noise r]{t), with noise strength Dr and correlation 
time Tr. In what follows we will describe the effect of 
external noise in terms of the effective parameters Dr 
and Tr in the rest of our analysis. 



III. THE FOKKER-PLANCK DESCRIPTION IN 
THE OVERDAMPED LIMIT 



In the over damped limit, Eq. reads 

i(t) : 



--V'ix) + -vit), 

7 7 



Wix)+^it), 



(22) 
(23) 



{nit)nit'))c = 



M!ieM'\it-t')\/Tc). (17) 



On the other hand, when the external noise correlation 
time is smaller than the internal one, we recover Eq. (|15p . 
In what follows, we shall focus on the situation when 
Tc » Tc- Thus, in terms of the effective noise ??(/:), the 
Langevin Eq.Q can be written as 



where we have defined Wix) = ~i\/^)V'ix) and ^(i) = 
(l/7)77(t). Clearly, ^(t) is the scaled noise and as -qit) is 
Gaussian [since fit) and e(t) are assumed to Gaussian] 
^(t) is also Gaussian. The Gaussian nature of e(t) is 
expressed by the probability distribution function 



P[C(i)] = 7Vexp 



--Ids 



ds'^is)as')pis - s') 



,(24) 



X = V, 

V = -V'ix)- I dt'-fit-t')vit') + r]it), (18) 
Jo 

which reduces to 

x = v,v^-V'ix)--fvit)+r]it), (19) 

where we have assumed that the internal noise /{t) is 
(5-correlated and the internal dissipation is Markovian so 



where P is the inverse of the correlation function of ^(t), 
and N is the normalization constant expressed by a path 
integral over ^(t). 



I I ds I ds'fis)fis')Pis - s') 



(25) 



Now, let i^it)) = and (C(t)^(t')> = cit~t'). Then from 
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Eq.dlSl) we get 

SN 



Now functional integration by parts yields 



+ I dsf{s)P{s - t) 

1 



ds'cit^s') J V^P[^]-^S{y - x{t)) 



^ I ds I ds'e(s)C(s')/3(s - s') 



X exp 

TV / dsP{t~s){^{s)) =0. 



(35) 



(26) 



Therefore, it follows that 

spm] 

Consequently, 
and 

= 



Again from Eq. ([23]) we have 



(27) 



ds J ds'fj{t - s)(3{t' - s')^{s)^{s') 



This equation possesses the unique solution 
dx(t) 

pm] 



e(t- i')exp 



dsW'{x{s)) 



(36) 



^2gyhere &{t — t') is defined by 



s'pim] 



smdmr 

ds J ds'l3{t - s)l3{t' - s')c{s - s') 
P{t-t'). (29) 



Q{t-t') = l,t>t' 
= 1/2, i = 
= 0,i < t'. 

Now substituting Eq.(l36l) into Eq.([35l) and from Eq. (|34l) 
we have 



Eq.dlll) imphes 

ds(3{t- s)c{s- s') = 5{t-s'), 



dP 
'dt 



(30) 



which shows that the kernel /3{s— s') is the inverse of the 
correlation function c(s — s'). Now, using Eqs. ((27|) and 
([30|) one may observe that 

Pimm) = -j ds'dt ~ s')^-^Ml. (31) 

The path integral is used to define the probability dis- 
tribution functional for x{t), then the solution of Eq, 
becomes 



X j V^P[i] exp 
xS{y-x{t))}. 



ds'c{t - s') 
dsW'{x{s)) 



(37) 



Piv.t)^ / V(,P[(]5{v-m)- 



From Eq.(132]) it follows that 



dP 
'dt 



--g:^S{y~ x{t))± 



(32) 



(33) 



Eq. ([37|) is not a Fokker-Planck equation. The second 
term cannot be reduced to a term containing P{y,t) be- 
cause of the non-Markovian dependence on x{s) for s < t. 
Fortunately, in our case c{t — t') ~ (ZJ^/rfl) exp(— |i — 
t'\/Tji), where Dr = Dr/^^, is an exponentially decay- 
ing function and for large 7 (as we are dealing with the 
over damped case) decays rapidly. We now change the 
variable t' = {t — s') and observe that 



^ ds'c{t- s')exp dsW'ix 



where x = dx/dt and can be replaced by the right hand 
side of Eq.jSSI) to get 

S = -I: / 'DiPmy-m)[w{x)+m], 



dt'c{t') exp 



dsW'{x{s)) 



— f dt' exp{-t' /ti^) 

TR. Jo 



dt 



dy 



X exp 



-^[w{y)p] -^J v^pmmy~x{t)) 



t'w'ixit))--itrw"{x{t))iit) 



(38) 



Neglecting the (i')^ term which can be shown to 
(34) be valid self consistently for small r/j under Markov 
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approximatiorif^^i^ we get 



/ dt'c{t')exp 
Jo 

/ at exp 

tr Jo 



Integrating Eq. (|42)) we have the expression of stationary 
probabihty distribution in terms of stationary current as 



dsW'{s) 



+t'W'{xit)) 



Pst{x) 



where 



P'R Jo 



Dr 



1 - TRW'{x{t)) ■ 



(39) 



""^"^^ DnJo G{x') 



(45) 



(46) 



for sufficiently large t. Substituting Eq.jSll) into Eq. fST]) 
we obtain the Fokker-Planck-Smoluchowski equation cor- 
responding to Eg. p^ as 



is the effective potential for the problem. Now applying 
the periodic boundary condition Eq. (|43p we have from 
Eq.gll) 



dP _ d r 1 

dt dx\^f 



1 



7 dx \l + {TRh)V"{x) 



P 



G{0)Pst{0) = 

^Jr 
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Jo 



2tt 



e^(^)da;. (47) 



where P = P{x,t), is the probability density of finding 
the particle at x at time t. Defining an auxiliary function 
G{x), Ea. (|in)l can be rewritten as 



(40]Using Eq. (|l^ and applying the normalization condition 
[Eq. ISD], we get 



dP 
'dt 



d_ r 1 

dx [■-f 



V'ix) + ^^G{x) 
7 ox 



P 



(41) 



where G{x) = + {tr/j)V"{x)]. 

IV. APPLICATION 

A. Solution under periodic boundary condition 

In this subsection we consider the dynamics of a Brow- 
nian particle moving in a periodic potential under a con- 
stant external force F. Then the above Fokker-Planck- 
Smoluchowski equation, Eq. (|41|1 reads as 



dP _ d r 1 

dt dx\^ 



V'{x)-F+^^G(x) 
7 ox 



P 



In the over damped limit the stationary current is given 
by 



J 



V'{x)~F+^^G{x) 
7 ax 



Pst{x), (42) 



where Pst{x) = P{x,t —>■ oo). Now under symmetric 
periodic potential with periodicity 27r, i.e., V{x) = V{x + 
27r), we may employ the periodic boundary condition and 
normalization over one period on Pgt (x) , 



P,tix)=P,tix + 2T:), 



(43) 



27r ^-Uix) 

'gW 



G(0)F,i(0)- r e'^^-'Ux' 
P'R Jo 



dx = 1. 



(48) 



Ehmination of G(0)Pst(0) from Eqs. (|47M8|) gives the ex- 
pression of stationary current 



J 



D, 



(l-e^(2.)) 



M 



(49) 



where 

M = 



dx I e^(^ - (l - 



G(x 

27r 







1 



Now, the average velocity, (v) — (i) , is given by 

(v) = ^((F-Vix))), 
7 

= - / {F~V'{x)}Pst{x)dx, 
7 Jo 

= - / 7^ + —TrG{x)Pst (x) dx, 
1 Jo { 1 ox J 

where we have made use of Eg.ip^. Since Pst{x) and 
V(x) are periodic functions of x, with period 27r, (w) is 
thus given by the constant probability current density, 
multiplied by 27r: (v) = 27r J. For the periodic potential 



2n \rt 1^1 



— — — — dx = -Fa say , 
G(x') 



and 



/ Pst{x)dx^l. 
Jo 



(44) 



with 



27r 



dx' 



DrJo G(x')' 
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By definition , the mobifity is given by 



1- ^ J 

II = um — — = 2n lim — . 

^ F^o F f^qF 



If we consider mobility only in linear response regime'^^, 
the double integral term in the expression of steady state 
current becomes 



J 



(1 -e--^«) 



M 



(50) 



where 

M : 



27r ^-U{x) 

dx 

G{x) JO 



eJJ^-'Ux' - (1 - e-^") 



2ir 



1 



~U{x) 



'a '^G{x) 

vanishes as F —i- and hence the mobility is given by 



(51) 



ir'^dxj,'-e^i^')dx' 
Using Einstein relation'^S, the diffusion rate is given by 



D 



47r2 
7 



(52) 



d4) 



teui^)dx!- 



G{x) 



-dx 



The above expression for diffusion rate is exact for any 
periodic potential and for any Gaussian noise process 
with decaying memory kernel. For a simple choice 
of potential, V{x), the above expression is analytically 
tractable. As for example, if we choose V{x) = Acosx, 
then to first order in t/j (assuming the damping is large) 



D = 



1 + Tfl(27-1)A 



^i(^) 



HA) 



(53) 



where 1^ is the modified Bessel function of order u. From 
Eg. ([55)1 ■ the diffusion rate is seen to increase from white 
to colored noise. The quantity A{Ii[A) / Io{A)) is strictly 
positive, except at A = where it is zero. The increase in 
D is consistent with the fact that the diffusion coefficient 
also increases from white to colored noise by a factor 
which is much larger than the enhancement factor of the 
potential. 



B. Bistable Kinetics 

The dynamics of a Brownian particle in a bistable po- 
tential models several physical phenomena^, and the 
standard form of the potential is given by. 



V{x) 



(54) 



which has symmetric minima at a; = -ii^Ja/b with an 
intervening local maxima at a; = of relative height Eb = 
a? /Ah. We consider that our system of interest is moving 
in a bistable potential of the form Eq. and is coupled 
to a bath which is modulated by an external Gaussian 
noise, e(<) with the statistical properties stated earlier. 
The corresponding Fokker-Planck equation in the over 
damped limit is given by Eq. (j40p . the explicit form of 
which with the potential as in Eq. ((54)) reads as 



dP_ 
'dt 




7 dx\^ + TR^{-ax^/2 + hxilA} 



P 



At steady state, the above equation reads as 



— i-ax + bx'^) Fix) 
dx ' 

. Dr d^ 



7 dx'^ 



7 



7 + TRi—a + 3bx^) 



Fix) = 0. (55) 



Since at steady state the stationary current vanishes. 
Eg. ([55]) takes the form 



where 



dF{x) 
dx 



R{x) 



+ R{x)F{x) ^ 0, 



g'{x) h{x) 



g{x) DRg{x) ' 



with 



7 - tr{o. - 36x2) 



, h{x) = ax — bx^ 



(56) 



(57) 



(58) 



and prime (') denotes differentiation with respect to x. 
The solution of Eq. ([TO]) is given by 



F{x) = N - TR{a - ibx^ 



exp 



2D, 



-(7 - aTR)x^ 



AabTR - 67 4 TRb g 
-X — „ „ x 



4L>f 



2Dr 



(59) 



where N is the normalization constant. From Eg. (j59p . 
we observe that the pre-exponential factor behaves like a 
constant in comparison to the exponential factor, which 
is an exact statement for tr = 0, i.e., when the external 
noise e(<) is (5-correlated. 

To this end, following the standard technique^^, the 
barrier crossing rate is obtained as 



y/2a 
7r7 



7 - o,tr 



7 + 2aTR 



exp 



0^7 
' AbDR 



2a 

7 



-TR 



,(60) 



which is valid in the strong friction regime. 
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FIG. 1; (color online) Plot of barrier crossing rate k as 
a function of dissipation constant 7 for various values of 
the effective correlation time tr. The solid lines are drawn 
from the theoretical expression, Ea. (|60p and the symbols are 
the results of numerical simulation of Eq. (|22|) . The val- 
ues of the parameters used are ksT — Dg = 0.05, Te = 
0.04 (red: solid), 0.06 (green: dash) and 0.08 (blue: dot) and 
Kq = 3 (red: solid), 5 (green: dash) and 7 (blue: dot). Inset: 
The normalized steady state PDF [ Ea. (f59|l ] using 7 = 5.0 
and the same parameter set as in the main figure. 



FIG. 2: (color online) Plot of barrier crossing rate k as 
a function of dissipation constant 7 for various values of 
coupling parameter k§. The solid lines are drawn from 
the theoretical expression, Ea. (l60p and the symbols are 
the results of numerical simulation of Ea. (|22p . The values 
of the parameters used are ksT = = 0.05 and Te = 
0.0625 (red: solid), 0.06 (green: dash) and 0.0583 (blue: dot) 
so that the effective correlation time tr always remains 0.05. 
Inset: The normalized steady state PDF [ Ea. (l59p ] using 
7 — 5.0 and the same parameter set as in the main figure. 



V. NUMERICAL IMPLEMENTATION 

To check the validity of our analytical result we nu- 
merically simulate the Langevin equation, Eq. (|22p us- 
ing Heun's algorithm^. In our simulation, we have al- 
ways used a small integration step At = 0.001 to ensure 
numerical stability. In addition to that all our numeri- 
cal results have been averaged over 10, 000 trajectories 
to obtain a smooth numerical profile. As mentioned in 
Sec. II, although we mention the explicit values of the 
external noise parameters (I?e and r^) used in the sim- 
ulation, we interpret our results in terms of the effective 
noise parameters Dfi and r^. 

In Fig.(l) we show the profile of escape rate, fc as a 
function of the dissipation constant, 7 for different val- 
ues of the effective correlation time, t^j. Numerically, 
the escape rate has been defined as the inverse of the 
mean first passage time^. The values of the different 
parameters used in the simulation are given in the figure 
caption. The profiles show that the numerical results are 
in good agreement with the analytical ones. As expected 
for a fixed value, the escape rate decreases with 7 but 
for a fixed value of 7, the escape rate increases with the 
increase of the effective correlation time, t^. To under- 
stand this behavior we numerically calculate the steady 
state PDF— for a fixed value of 7, which is a measure of 
the dynamics in the bistable potential. The analytical [ 
Eq. ([55)1 ] and numerical profiles of steady state PDF [see 
the inset of Fig. (1)] show that the barrier height decreases 
as Tji increases which effectively increases the escape rate, 
k. In Fig. (2) we show the variation of the escape rate 
with different values of the coupling term kq which re- 
fiects an increasing trend of k with kq. The steady state 



PDF profile accounts for this behavior with decrease in 
the barrier height [see the inset]. 



VI. CONCLUSION 

A system reservoir model, where the reservoir is modu- 
lated externally by a Gaussian colored noise has been pro- 
posed to study the transport of an overdamped Brown- 
ian particle in a periodic potential. Based on the Fokker- 
Planck-Smoluchowski description we calculate the mobil- 
ity of the Brownian particle in the linear response regime 
and using Einstein's relation the diffusion rate is calcu- 
lated for any arbitrary periodic potential where the ex- 
ternal driving noise is colored. For a cosine potential we 
obtain the diffusion rate in a closed analytical form and 
observe that the diffusion rate increases from white to 
colored noise. As an immediate application of our for- 
malism, we study the bistable kinetics of the Brownian 
particle and demonstrate the dependence of the corre- 
lation time of the external colored noise, by which the 
bath is modulated, on the steady state probability den- 
sity function and observe the barrier crossing dynamics 
to obtain the expression for escape rate. Our analyti- 
cal result for the escape rate is then compared with the 
Langevin simulation result which shows that both are in 
very good agreement. To the end it should be noted that 
nonlinear system reservoir coupling may also be consid- 
ered to obtain a Langevin equation for a Brownian parti- 
cle effectively driven by a state dependent colored noise, 
from which one may observe various dynamical and kine- 
matical aspects of the Brownian particle. We hope to 
address these issues in near the future. 
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